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ABSTRACT 

Accurate and radially extended stellar kinematic data reaching R = 97" from the 
center are presented for the cD galaxy of Fornax, NGC 1399. The stellar rotation 
is small (< 30 km/s); the stellar velocity dispersion remains constant at 250-270 km/s. 
The deviations from Gaussian line of sight velocity distributions are small, at the percent 
level. 

We construct dynamical models of the galaxy, deprojecting its nearly round (EO-El) 
surface brightness distribution, and determining the spherical distribution function that 
best fits (at the 4 percent level) the kinematic data on a grid of parametrized potentials. 
We find that the stellar orbital structure is moderately radial, with f3 = 0.3 ±0.1 for 
R < 60", similar to results found for some normal giant ellipticals. The gravitational 
potential is dominated by the luminous component out to the last data point, with a 
mass-to-light ratio M/Lb = 10 Mq/Lq, although the presence of a central black hole 
of M ^ 5 X IO^Mq is compatible with the data in the inner 5 arcsec. The infludence 
of the dark component is marginally detected starting from R ^ 60". Using the radial 
velocities of the globular clusters and planetary nebulae of the galaxy we constrain the 
potential more strongly, ruling out the self-consistent case and finding that the best fit 
solution agrees with X-ray determinations. The resulting total mass and mass-to-light 
ratio are M = 1.2 - 2.5 x lO^^M© and M/Lb = 22 -48 Mq/Lq inside R = 417" or 35 
kpc for D=17.6 Mpc. 
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and lenticular, cD - galaxies: individual: NGC 1399 - galaxies: kinematics and dynamics 
- dark matter. 
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1. Introduction 

Determining the gravitational potential of ellipti- 
cal galaxies and the amount of dark matter present 
is a notoriously difficult problem. This stems from 
the absence of extended cold gas which could be used 
as dynamical tracer, and the unknown intrinsic shape 
and orbital structure (e.g., Saglia 1996 and references 
therein). Progress has been made in recent years 
thanks to improved observational data and dynam- 
ical modeling techniques. 

On the one hand, we are now able to measure not 
only the stellar velocity and velocity dispersion pro- 
files, but also the shape parameters of the line of sight 
velocity distributions with good accuracy (van der 
Marel and Franx 1993; Bender, Saglia and Gerhard 
1994, hereafter BSG94; CaroUo et al. 1995; Gerhard 
et al. 1998a, hereafter G-f98; Statler and Smecker- 
Hane 1999). Such data contain information on the 
anisotropy and the gravitational potential (Gerhard 
1993, Merritt 1993). On the other hand, algorithms 
based on Schwarzschild's (1979) orbit superposition 
technique (Rix et al. 1997; Cretton et al. 1999) 
or on the non-parametric reconstruction of the un- 
derlying distribution function ( G-f 98; Matthias and 
Gerhard 1999) allow a detailed modeling of the data. 

Up to now, only a limited number of ellipticals have 
been analyzed with these methods (Rix et al. 1997; 
G+98; Matthias and Gerhard 1999). Nevertheless, 
the following pattern appears to emerge from these 
studies. Out to 1-2 half-luminosity radii Re these 
objects are moderately radially anisotropic, and they 
have a dark component which starts to dominate the 
dynamics at about R > 1 — 2_Re. This is consis- 
tent with previous studies of extended velocity dis- 
persion profiles of ellipticals (Saglia, Bertin and Sti- 
avelli 1992; Saglia et al. 1993), of the radial velocities 
of globular clusters (Mould et al. 1990; Grillmair 
et al. 1994; Cohen and Ryzhov 1997; Minniti et 
al. 1998, Kissler-Patig et al. 1998) and planetary 
nebulae (Arnaboldi et al. 1994, 1998; Hui et al. 
1995); of the X-ray emission around ellipticals (Ikebe 
et al. 1996, Jones et al. 1997) and, more recently, 
the modeling of strong lensing systems (Romanowsky 
and Kochancck 1999) and the statistics of galaxy- 
galaxy lensing (Geiger and Schneider 1998). 

In this paper we apply the non-parametric distri- 
bution function reconstruction method to new, accu- 
rate and extended data of the Fornax cD galaxy NGC 
1399, to explore whether the conclusions reached for 



normal ellipticals are valid also for the class of the 
most massive early- type galaxies. A preliminary de- 
scription of the results can be found in Kronawitter 
et al. (1999). In the following we adopt a distance 
of 17.6 Mpc. NGC 1399 has been the subject of a 
number of dynamical studies in the past. Bicknell 
et al. (1989) analyzed the major axis velocity dis- 
persion profile using the Jeans equations, concluding 
that a constant mass-to-light ratio model with little 
anisotropy is able to reproduce the data out to ~ 85". 
More recently, Graham et al. (1998) reached the 
same conclusion using models based on i?^/" lumi- 
nosity profiles. X-ray observations indicate that the 
outer parts (i? > 120" = 10.2 kpc) are dominated by 
a dark halo which merges into the one of the Fornax 
cluster at around i? « 53 kpc (Ikebe et al. 1996). 
However, the potential inside i? « 10 kpc is little con- 
strained due to the presence of a cooling flow (Ran- 
garajan et al. 1995). The radial velocities of globular 
clusters and planetary nebulae (see references above) 
agree with these mass determinations. 

The structure of the paper is as follows. We present 
the data in §H, where also the data reduction and the 
comparison with the literature are described and dis- 
cussed. We review our dynamical modeling approach 
in §H and discuss some modifications needed for the 
analysis of NGC 1399. Results are presented in 
where the constraints from the absorption line data 
are combined with those from globular cluster and 
planetary nebulae radial velocities and X-ray data. 
Conclusions are drawn in 

2. Observations of NGC 1399 

The observations were performed from the 6th to 
the 8th of December 1994, using EMMI and the NTT, 
in remote observing mode from Garching. The grat- 
ing # 6 (13.5 A/ mm) was used in combination with 
the ESO CCD Tek #36 (2048^ 24 //m pixels, 0.268 
arcsec/pixel) and a 3 arcsec wide slit, giving a FWHM 
spectroscopic resolution of 3.5 A. The central wave- 
length was 5150 A, giving 85 km/s (Jinstr resolution; 
the spectral range (4826-5474 A) covered the absorp- 
tion lines H/3, Mg and Fe. A 1.5 hour long spec- 
trum was taken along the major axis of the galaxy 
(PA=110°); additional 9 hours of integrations were 
taken parallel to the minor axis, shifting the slit 42 
arcsec from the center. Numerous template stars were 
also observed, wiggled and trailed along the slit. 

The standard reduction steps (bias subtraction, 
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flat fielding, cosmic ray removal, logarithmic wave- 
length calibration, sky subtraction) were performed 
under MIDAS. A sky subtraction better than 1 per 
cent was achieved. The analysis of the data was car- 
ried out using the FCQ method (Bender 1990) follow- 
ing BSG94 and G-l-98 . The spectra were rebinned 
along the spatial direction to obtain a nearly constant 
signal-to-noise ratio larger than 50 per resolution ele- 
ment. The effects of the continuum fitting and instru- 
mental resolution were extensively tested by Monte 
Carlo simulations. The residual systematic effects on 
the values of the and /i4 profiles are less than 0.01, 
and less than 1.5% in a. The errorbars (which reflect 
the random errors and do not take into account sys- 
tematic effects such as template mismatching or the 
presence of dust and faint emission, see below) are 
very small, in the range 3-7 km/s for the recessional 
velocities and velocity dispersion, and 0.006-0.02 for 
the /13 and /14 coefficients. They are calibrated with 
simulations and are determined to better than 20%. 
They agree with the scatter of the points except for 
the 5 pairs of points with 5 < i? < 20 arcsec. In this 
region the spectra on the two sides of the galaxy are 
slightly different, probably due to low level emission. 
The Ha maps (Goudfrooij et al. 1994, Singh et al. 
1995, Macchetto et al. 1996) show that the emission 
is patchy in the 10-20 arcsec region along the slit po- 
sition. Dust is also present (Goudfrooij et al. 1994). 
Bicknell et al. (1989), who notice the same asym- 
metry, suggest that "old shells" might also be partly 
responsible. In the inner 1.5 arcsec, where a increases 
to 370 km/s, there is still some template mismatch- 
ing, which none of the 15 comparison stars is able to 
eliminate. The outermost datapoint has low signal 
and only the mean velocity dispersion datapoints are 
reliable. 

Fig. shows the kinematical profiles including 
and /i4, with the major axis data folded with respect 
to the center, and the spectra parallel to the minor 
axis folded with respect to the major axis. The rota- 
tion along the major axis is small, <J 30 km/s, and is 
present also along the position parallel to the minor 
axis. The velocities measured in the inner 3 arcsec 
might indicate the presence of a kinematically decou- 
pled core. The velocity dispersion increases to « 370 
km/s in the inner 5 arcsec, and fiattens to w 250 
km/s in the outer region. The antisymmetric devia- 
tions from a Gaussian profile as parametrized by the 
/13 coefficients are slightly negative where rotation is 
detected (see BSG94 ). The symmetric deviations 



measured by the /14 coefficients are smaller than 0.05 
and mostly positive. 

The comparison with the rotational velocities along 
the major axis of Bicknell et al. 1989), D'Onofrio et 
al. (1995), Franx et al. (1989), Graham et al. (1998) 
and Longo et al. (1994) is shown in the lower panel 
of Fig. ||. There is overall good agreement within the 
quoted errors. The systematic deviations observed for 
the outermost datapoints of Longo et al. (1994) and 
D'Onofrio et al. (1995) suggest that these datasets 
are probably affected by residual sky subtraction er- 
rors. The comparison with the dispersion profiles of 
Bicknefl et al. (1989), D'Onofrio et al. (1995), Franx 
et al. (1989), Graham et al. (1998), Longo et al. 
(1994), StiaveUi et al. (1993), and Winsall and Free- 
man (1993) is shown in the upper panel of Fig. |[ 
There is overall agreement within the (usually large) 
errorbars given by the authors. Systematic differences 
are observed in the inner 5 arcsec, where differences 
in seeing and slit width may play an important role. 
In the outer parts the velocity dispersions of Bick- 
nell et al. (1989) are systematically smaller than our 
data. However, our data agree well with Graham et 
al. (1998) and with the point at i? w 80 arcsec mea- 
sured by Winsall and Freeman (1993). 

3. Dynamical modeling 



The algorithm used to analyse the kinematic data 
presented above follows closely the technique dis- 
cussed in detail by G-l-98. Therefore here we limit 
ourselves to summarizing its main features and the 
modifications we implemented to fit NGC 1399. In 
§3.1 we f ocus on the modeling of the light distribu- 
tion; in §3.2 we concentrate on the non-parametric 
reconstruction of the distribution function. 



3.1. The deprojection of the light distribu- 
tion 

The surface-brightness (SB) data of NGC 1399 
consist of three sets: in the innermost part we use 
V-band HST-data from Lauer et al. (1995), the 
medium range is covered by ground-based B-band 
CCD-data from Bicknell et al. (1989), and the outer 
part by photographic V-band data from Schombert 
(1986). As described in Bicknell et al. (1989) we used 
the Schombert (1986) data shifted by 1.2 mag and 
from 57" on outwards. The HST-data are also shifted 
by 0.8 mag and are used out to 10". With an effec- 
tive radius of 42" according to Faber et al. (1989) the 
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Fig. 1.— The kinematics of NGC 1399. (a) The folded 
mean velocity, (b) velocity dispersion, (c) ha and (d) /i4 
profiles. Crosses and filled circles refer to the two sides of 
the galaxy and the major axis spectrum. Open and filled 
squares refer to the two sides of the galaxy and the spectra 
taken parallel to the minor axis and shifted 42 arcsec from 
the center. 




Fig. 2. — Comparison with the velocity dispersion (top) 
and rotation velocities (bottom ) of Bicknell et al. (1989, 
crosses), Franx et al. (1989, triangular diamonds), Gra- 
ham et al. (1998, open squares), Longo et al. (1994, open 
triangles), D'Onofrio et al. (1995, triangular crosses), Sti- 
avelli et al. (1993, stars), Winsall and Freeman (1993, 
open circle). The line connects the data points determined 
here (full circles). Note the change of scale at i? = 20". 
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photometric data extend over 75 effective radii. (Note 
that Caon et al. 1994 find = ^/aJTe = 127" by fit- 
ting their extended photometry. The determination 
of the total magnitude and effective radius of NGC 
1399 is highly uncertain, since its hght profile is very 
shallow in the outer regions. The exact value of Re 
does not play a role in the following analysis). 

The surface brightness profile deviates strongly 
from a de Vaucouleurs or Jaffe profile, being much 
brighter in the outer parts. Thus contrary to NGC 
6703, for which G+98 could model the light distri- 
bution of NGC 6703 as a Jaffe (1983) profile, we 
must here determine the three-dimensional luminous 
density of NGC 1399 by deprojection. 

NGC 1399 is nearly round (EO-El). If we assume 
spherical symmetry the deprojection of the SB profile 
E(i?) to a spatial luminosity density ^'""(7-) is unique 
and is given by the standard Abel equation 
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To avoid the amplification of noise, we follow Wahba 
and Wendelberger (1980), Scott (1990), and Geb- 
hardt et al. (1996), and first smooth the data by find- 
ing the function S(-R) (expressed as thin plate splines) 
which maximizes the logarithm M of likelihood C 



M = \ogC 



(logs,- log I](i?,))' 



AP(E), (2) 



where Si denotes the observed SB at projected radius 
Ri and the measurement error. The second term on 
the right hand side is a penalty function P (the sec- 
ond derivative of the SB profile — 2 . log S/d log R^ ) 
which penalizes oscillations in S(7?i) induced by noise. 
The factor A controls the degree of smoothing. We 
fixed it to 0.003, which allows to reproduce smoothly 
the outer, noisy regions of the galaxy without affect- 
ing the inner, high signal to noise datapoints. We also 
tried using the more objective Generalized Cross Val- 
idation (GCV) method to determine the smoothing 
parameter, but this did not give useful results for this 
inhomogeneous data set. 

The smoothed SB-profile is then deprojected using 
eq. (1). The error introduced by this inversion can be 
checked by generating data from a known model and 
deprojecting. For SB data drawn from a Jaffe-model 
(Jaffe 1983), with similar quahty as the NGC 1399 
data, the deviations between the known and the de- 
projected luminosity density profile are smaller than 



0.5% rms. Figure || shows the surface brightness pro- 
file of NGC 1399 and that obtained by projecting the 
luminosity density model computed from our depro- 
jection method. The differences are smaller than 0.05 
mag/arcsec^ for most of the radial range, and smaller 
than 0.02 mag in the integrated luminosities. 

The resulting density-profile is rather shallow, with 
some intrinsic structure at large radii. This would 
be lost if we had fitted a parametrized model to the 
surface brightness data; by contrast the regression- 
inversion method described above does not smooth 
away intrinsic features or force the profile to a fixed 
asymptotic decrease in the outer parts, as is the case 
in many analytic models. 

3.2. The distribution function 

To analyse the line of sight velocity profile data 
we proceed in two steps (Merritt 1993, Rix et al. 
1997, G-l-98 ): First, the gravitational potential is 
kept fixed and the best-estimate distribution function 
(DF) is determined, in a least -squares sense. For ideal 
data, the solution to this problem in the spherical case 
would be unique (Dejonghe and Merritt 1992). Then, 
this process is repeated for different gravitational po- 
tentials to determine those for which a statistically 
good fit to the data can be achieved with positive 
DF. Because changing the gravitational potential has 
a different effect on the line profiles than changing 
the anisotropy of the stellar orbits (e.g., Gerhard et 
al. 1998b), both can be constrained; i.e., only a lim- 
ited range of potentials is consistent with the data. 

The potentials examined here are the sum of the 
potential of the luminous matter, obtained by solving 
Poisson's equation for the spherical luminosity den- 



sity found in §3.1, and that of a dark halo taken to 
be 



DM 
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(3) 



where is the asymptotic circular velocity and ro 
the core radius of the halo. We explored models with 
velocities wq from (the self-consistent model) to 780 
km/s and ro from 0.5 to 6.0 Re- Models with these 
potentials minimize the contribution of the dark halo 
near the center and hence maximize the mass-to-light 
ratio of the luminous component. In this sense they 
are similar to the maximum disk models for spiral 
galaxies; by analogy we refer to them as maximum 
stellar mass models. A priori these halo potentials 
might not be correct; in fact, from cosmological sim- 
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Fig. 3. — Top panel: the combined surface brightness 
profile fi of NGC 1399 (crosses) and the projection of the 
derived luminous density model (full line; see §3.1). Mid- 
dle panel: the residuals A/i in magnitudes per square arc- 
sec between the observed SB profile and the model. Lower 
panel: the residuals Amag in magnitudes between the in- 
tegrated luminosity of NGC 1399 and the model. 



ulations (Navarro, Frenk and White 1997) one ex- 
pects more centrally peaked halo distributions with 
asymptotic profiles at large radii oc r""^. However, the 
mass profiles derived from our approach match well 
the X-ray mass distribution (see Fig. ^ and §4.3). In 
addition, the luminous density profile derived in §3.1 
is well approximated by a r^^ power law beyond 1" . 

To recover the DF of a galaxy we construct a set 
of basis DFs in a fixed potential. Each basis DF has 
a different anisotropy structure and reproduces the 
luminous density profile determined as in §3.1. The 



total DF is then a weighted sum of 57 of these basis 
functions: 



(4) 



with flfc = 1 and the additional constraint that 
J2k O'kfk > 0. everywhere. For the basis functions fk 
we used sequences of tangentially anisotropic models 
as described in §4.1 and App. A of G-l-98 , as well 
as the isotropic model and a function with radially 
anisotropic orbital structure. The radial function is 
used for determining the optimal velocity scales when 
computing the projected Gauss-Hermite moments for 
each of the fk- The coefficients are determined 
by fitting the projected model kinematics to the / 
datapoints such that 



E 

k=l,I 



is minimal. Here xi-i xi ^^'^ xl ^re as in Eqs. (8) 
and (9) of G+98 and measure the error-weighted 
square differences between the velocity dispersions of 
the galaxy and the model, and similarly for the sec- 
ond and fourth order Gauss-Hermite moments related 
to the measured hi coefficients. Following G-l-98 , we 
write the penalty function A as a combination of the 
second derivatives of the composite DF, evaluated on 
a grid of points {Ei,Xj) in the energy and circular- 
ity integrals. However, in the present case some care 
is needed in normalizing these second derivatives, to 
avoid that a few data points dominate the regularisa- 
tion terms. For, as expected from the correlation be- 
tween galaxy magnitude and the steepness of the cen- 
tral profile slope (Faber et al. 1997), NGC 1399 has 
a rather shallow central cusp. Flat density-profiles 
towards the center cause the distribution function to 
become very steep (Dehnen 1993), and similarly its 
second derivatives. After testing different regulariza- 
tion functionals on simulated data of kwown distri- 
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bution functions, we found that the best results are 
obtained by setting 

I dEdx ) 



A(/) 



1 



^2 \dE^ ) 



1 



where /iso, /4o //so ^'^^ isotropic distribution func- 
tion and its first and second derivative with respect 
to energy. 

The value of A in Eq. ^ determines the degree of 
smoothness of the composite DF. Low values of A re- 
sult in a perfect fit to the data at the price of unphys- 
ical small-scale variations in the DF. Large values of 
A produce very smooth DFs which might not give ac- 
ceptable fits to the data. We estimate the optimal 
value of A by considering the composite DF recovered 
from Monte Carlo simulated data, which are derived 
from an underlying distribution function similar to 
the one that gives the best fit to the NGC 1399 data, 
with the assumption of Gaussian errors (see G-l-98 
for details). Fig. ^ shows the normalized total of 
the fit to the artificial data and the percentage rms 
variations of the recovered DF, averaged over 25 re- 
alizations, as a function of A. A potential with only 
the luminous matter and one including a halo with 
To = 168 arcsec, vq — 427 km/s are used for Fig. ^. 
In both cases values around A w 0.01 give « 1 
and less than 15% mean rms variations of the DF. In 
the analysis of NGC 1399 we thus adopt A = 0.02. 
Note that in G-l-98 we found that slightly different 
values of A were needed to fit the self-consistent and 
the halo potential cases. The difference here stems 
from the dominant role of the luminous component 
in NGC 1399 (see 

Finally, we constructed cumulative distribu- 
tions for sets of 100 Monte Carlo simulations of our 
data set, drawn from radial and tangential models 
in self-consistent or halo potentials, respectively, and 
analyzed using A = 0.02. We find that the 95% confi- 
dence level is reached at a value of ~ 1-3 per data 
point in all these cases. The same value results also 
when all errors are increased by 25%, with A = 0.02 
still beeing the optimal smoothing parameter. If we 
increase all errors by 40%, the 95% confidence level is 
reached at a value of « 1.37 per data point. 
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Fig. 4. — Model results as a function of the regularization 
parameter A. The upper two curves show the total per 
o and hi data point of regularized composite models with 
57 basis functions. The lower curves show the rms devia- 
tion between the recovered composite model distribution 
function and the true distribution function from which the 
data sets were generated. This rms deviation was evalu- 
ated on a grid in energy and angular momentum corre- 
sponding to radii < 3 times the radius of the outermost 
data point. Triangles refer to the self-consistent model, 
pentagons to a model in a halo potential with ro = 168 
arcsec and vo = 427 km/s. 



4. Results 

4.1. Stellar kinematics 



As discussed in §|[ the kinematic data show differ- 
ences between the two sides of the galaxy significantly 
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larger than the (small) statistical errors, probably due 
to low-level patchy emission. For radii 10 < |i?| < 22 
arcsec, we therefore use in the following the average 
values of a and over the two sides, adopting as 
errors half of the side-to-side differences. In addi- 
tion, while minimizing the square differences between 
models and datapoints over the whole radial range, we 
quote values computed only from points at \R\ > 5 
arcsec. The velocity dispersion of NGC 1399 steeply 
increases to 370 km/s in the central 5 arcsec, which 
probably indicates the presence of a central massive 
object (Stiavelli et al. 1993). We explore models 
with central black holes at the end of this section. 

The full line in Fig. |^ shows a dynamical model 
that best fits the kinematic data, obtained with a dark 
halo potential with ro = 210 arcsec, vq — 323 km/s, 
corresponding to a total circular velocity Vdvc — 383 
km/s at the last data points, R — 97". Fits of sim- 
ilar quality are obtained for a small range of model 
potentials; see the discussion in G-l-98. Fig. || shows 
three other models that span this range in the present 
case. These have circular velocities at the last data 
point of Vcirc = 404,425,442 km/s, respectively. All 
these models represent the data very well, with rms 
differences in a of about 10 km/s, or 4.4 %, and in 
/i4 of about 0.018 (excluding the last datapoint). The 
dotted line in the figure shows the best model with- 
out dark halo {Vcn-c — 364 km/s); this has similar rms 
differences in a and when only datapoints within 
60 arcsec are considered, but substantially larger val- 
ues, 19 km/s, or 8.7% rms in a, and 0.027 rms in 
/i4, when the datapoints at i? > 60" are included. 
Therefore we conclude that the gravitational poten- 
tial of the galaxy in the region probed by the absorp- 
tion line data is dominated by the luminous compo- 
nent, but that the influence of the dark component is 
marginally detected starting at i? « 60". 

Despite the excellent fit, the statistical interpreta- 
tion of the resulting values is not straightforward. 
The top panel of Fig. || shows the values of the nor- 
malized global xi+hi obtained for the model fits, as 
a function of the circular velocity of the model at the 
last datapoint, R — 97". Best-fit models are found 
in the range 380 < Vcirc < 450 km/s with values 
xi+hi ~ 1-9' while the self-consistent model gives 
xi+hi ~ 3.7. The Monte Carlo simulations of our 
data shown in Fig. ^ demonstrate that the employed 
set of basis functions is sufficient to produce « 1, 
for a Gaussian error distribution and when the under- 
lying model distribution function is similar to NGC 



1399 and smooth. In these simulations the 95% con- 
fidence limit corresponds to ~ The rather 
high values obtained for the real dataset might 
partly be caused by an underestimation of the error 
bars. Fitting to the data for NGC 1399 with all error- 
bars increased by 25% gives a normalized = 1-3 for 
our best-fitting model, just marginally acceptable in a 
statistical sense (see above) . Only when all errorbars 
are increased by 40% is the fitted model's reduced 
to 1.06. However, the (statistical) uncertainty of the 
Monte Carlo-estimated error bars is not larger than 
20% (see §^ , so most of the large x^ must be due to 
systematic effects. 

Smaller, but still greater than unity, values of x^ 
are obtained if the two sides of the galaxies are fit 
separately (x^ = 1-19 and 1.58 for the two sides, us- 
ing the original errors). This indicates that the sys- 
tematic differences between the two sides discussed 
in §1^ are playing an important role. Also important 
are sudden point-to-point variations in the measured 
data values which are unlikely to be physical. We 
conclude that the residual differences (at the two- 
three percent level) between model and data are of a 
systematic nature, with the most obvious candidates 
being dust, unresolved emission, local template mis- 
match, and slight oversmoothing of the model. This 
also suggests that we might have reached the intrin- 
sic precision limit of stellar absorption line measure- 
ments. 

To estimate confidence intervals for the potential 
parameters, we therefore proceed as follows. We as- 
sume that due to systematic effects the real errors on 
our kinematic data points are 40% larger than the 
statistical errors determined from Monte Carlo simu- 
lations, so that the best-fitting model has a x^ — 1 
per data point (the number of degrees of freedom is 
only slightly smaller than the number of data points 
due to the enforced regularization) . The 95% confi- 
dence line is then at x^ = 1-37 x 1.4^ = 2.68 in units 
of the original errors; see Figure ^ top panel. The re- 
sulting 95% confidence range for the circular velocity 
at i? = 97" is Fcirc = 420 ± 40 km/s. This is consis- 
tent with the sharp increase of the minimum values 
of x^ seen for Vdic < 380 km/s and Kirc > 450 km/s, 
which clearly indicates that this procedure gives a fair 
estimate for the range of allowed circular velocities for 
the combined luminous matter and halo potential. In 
this sense the four halo models shown in Fig. || bracket 
the range of acceptable models. 

The bottom panel of Fig. ^ shows the anisotropy 
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profiles derived for the different potentials discussed 
above. As a general trend, the profiles favour values 
oi (3 ~ 0.3. The maximum at 10" and sharp drop 
in the inner 5 arcsec is uncertain (see discussion be- 
low), while the tendency of the most massive models 
to become tangentially anisotropic in the outer parts 
is also uncertain until confirmed by data at still larger 
radii; see §4.2 and the discussion in G-l-98 . We con- 
clude that NGC 1399 is slightly radially anisotropic 
with (3 — 0.3 ± 0.1 out to i? = 60 arcsec or 5 kpc; 
and that the anisotropy is not yet well-constrained 
at larger radii. This is similar to what was found for 
NGC 6703 (G+98 ) and NGC 2434 (Rix et al. 1997). 

Finally, we briefly investigate the influence of the 
possible presence of a massive black hole on the in- 
ferred anisotropy in the center, without attempting 
a detailed analysis (this will be worthwhile only for 
data with HST resolution). We construct models with 
a central black hole of mass around M — 5 x lO^M©, 
which is approximately the mass expected for a galaxy 
of the luminosity of NGC 1399 (van der Marel 1999). 
We find that in these models the fit to the velocity dis- 
persion profile in the inner 5 arcsec follows the steep 
increase of the measured a better, and the anisotropy 
change is milder. In models with massive black holes 
the inferred anisotropy can be substantially changed 
even at 10", and is slightly reduced at ^ 20". 

4.2. The Globular Cluster and Planetary Neb- 
ulae systems 

The mass distribution of NGC 1399 can be con- 
strained further by considering its globular cluster 
(GC) and planetary nebulae (PN) systems. Grillmair 
et al. (1994), Minniti et al. (1998) and Kissler- 
Patig et al. (1998) measured the radial velocities of 
74 GCs, with typical errors of 100 km/s. Forbes et al. 
(1998) studied the radial distribution of the GCs and 
found that their surface density profile is fiatter than 
the galaxy starlight, following a power law of slope 
7 « —1.47, 0.28 dex shallower than the galaxy. Thus 
at the same radius we expect the GCs to have a veloc- 
ity dispersion ^ (1.75/1.47)^/^ « 1.08 larger than the 
stars of NGC 1399 (see, e.g., Binney and Tremaine 
1987), assuming a logarithmic potential and a nearly 
isotropic distribution function for both components, 
which for the stars is approximately the case. The dis- 
cussion below does not change if we ignore this small 
correction. 

Arnaboldi et al. (1994) measured the radial veloc- 
ities of 37 PNs, with typical errors of 50 km/s. The 
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1 10 100 

R [arcsec] 

Fig. 5. — Dynamical models for the kinematics of NGC 
1399 in several luminous plus dark matter potentials, com- 
pared to the measured projected velocity dispersion (top 
panel) and velocity profile shape parameter (middle 
panel). Full, dashed, long-dashed, dot-dashed lines: best- 
fit models with Vcirc = 383, 404, 425, 442 km/s at R = 97" 
respectively. Dotted lines: self-consistent model with lu- 
minous matter only, Kirc = 364 km/s. Filled squares re- 
fer to the stellar kinematics, open squares are the velocity 
dispersions deriv ed f or globular clusters, reduced by the 
factor 1.08 (see §4.2). The open triangles are the veloc- 
ity dispersions derived from planetary nebulae velocities. 
The bottom panel shows the models' intrinsic anisotropy 
parameter P{r). 
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intrinsic density distribution of the PNs is not as well 
established as for the GCs. We assume that in the 
outer parts of NGC 1399 the density distribution of 
PNs follows that of the stars, with the same veloc- 
ity dispersion. This is consistent with the histogram 
shown in McMillian, CiarduUo and Jacoby (1993), 
but the number of PNs used is small and there is a 
possible selection bias from the coupled effects of the 
metallicity dependence of the specific density of PNs 
(CiarduUo and Jacoby 1992) and the observed color 
gradients in elliptical galaxies (Richer, McCall, Ari- 
moto 1997). Fig. || shows the velocity dispersions 
(with their statistical errors) derived by binning radi- 
ally the GC and PN velocities. 

We proceed to compute the likelihood of the two 
datasets on our grid of gravitational potentials as fol- 
lows. For each model we take the projected velocity 
dispersion profile <j{R), as determined in §4J, and ap- 
proximate the probability for measuring a given line- 
of-sight velocity at radius Ri by a Gaussian of zero 
mean and rnis= fa{Ri), where / = 1.08 for the com- 
parison with the GCs and / = 1 for the PNs. This 
simplifying assumption is justified by the low number 
of GCs and PNs velocities available and their rather 
large errors. Monte Carlo simulations showed that 
deviations from Gaussians can be detected only with 
samples of more than a few hundred measured veloc- 
ities; see also Merritt (1997). After convolution with 
the (Gaussian) error distribution specified by the re- 
spective fJi of each data point, we write the likelihood 
of the set of measured radial velocities Vi as 



C 



GC 



n 



74 



1 



^2T:[pa^{Ri) + ai 



: exp 



2[/V2(i?0+^?] 
(7) 



and similarly for £pN- For comparison we also com- 
pute the likelihoods £mc of simulated sets of GC 
and PN radial velocities, generated from our best-fit 
model of NGC 1399. 

The middle and bottom panels of Fig. ^ show the 
results. Here we plot as a function of the circular 
velocity at i? = 97" the likelihoods of the models. 
Additionally, we show two lines. The dotted line dis- 
plays the mean value of the likelihood — ln£ for 100 
Monte Carlo realizations of iV = 74 GCs with veloc- 
ities drawn from the projected kinematics of one of 
our best fitting models (-547.6), and the same quan- 
tity (-266.45) for iV = 37 PNs. The dashed hue 
gives the 95% confidence level below which only 5% 
of the Monte Carlo realizations fall, for the GC (- 



553.9) and the PN (-272.5) samples. It is seen that 
the GCs and PNs favour the high mass halo models of 
the range compatible with the absorption line data, 
and are inconsistent with the self-consistent model. 
The models that are best for the GC sample have 
— In^GC ~ —547; they therefore give a good repre- 
sentation of the available data set. The same is true 
for the PN set. For both the GCs and the PNs, Cpn 
decreases sharply for models with T4irc < 380 km/s, in 
agreement with the results of §4.1. The Monte Carlo 



simulations quoted above rule out the self-consistent 
model to almost certainty with either the GC or the 
PN sample. However, models with more massive ha- 
los than allowed by the stellar kinematics are com- 
patible with the GC and PN data. 

We have estimated the effect of a possible contam- 
ination of the GC and PN samples by intracluster 
objects. As suggested by Theuns and Warren (1997) 
and Mendez et al. (1997), some of the GCs and 
PNs considered above could in fact be foreground or 
background objects in the Fornax Cluster but not as- 
sociated with NGC 1399. To test the sensitivity of 
our results to such a contamination, we repeat the 
analysis on a subset of GCs and PNs where the 10 
or 20 percent objects with the highest relative veloc- 
ities have been eliminated. This is an extreme case, 
since we might expect some of the contaminors to 
have smaller velocities. We find that the results de- 
rived for the GC set arc rather robust, with the best- 
fit potential remaining within the error range given 
above. However, when considering the PNs without 
the 8 fastest objects dark halos with circular veloci- 
ties as low as Vdrc ~ 350 km/s are allowed. Velocities 
Vcirc > 550 km/s are also less favoured. 

4.3. The mass and M/L distribution 

Fig. shows the permitted range for the cumula- 
tive mass distribution in NGC 1399, as derived from 
our combined fits to the absorption line kinematics 
and the radial velocities of GCs and PNs. The dark 
matter contribution is negligible for i? < 60 arcsec or 
5 kpc and becomes progressively dominant at larger 
R. It comprises at most 1/3 of the total mass at 
the last stellar kinematic data point, R — 97" or 8 
kpc, and 1-3 times the luminous mass at the radius 
of the most distant GC with a measured radial veloc- 
ity, R « 417 arcsec or 35 kpc. Inside this radius, the 
allowed mass range matches the profiles derived from 
X-ray data by modeling the cooling fiow (Ikebe et al. 
1996), giving M = 1.2-2.5 x lO^^AfQ. At larger radn 
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Fig. 6. — Quality of tlie kinematic model fits for NGC 
1399 as a function of the true circular velocity at _R = 97" . 
Top: the normalized of the fit to the a — hi datapoints. 
The dotted line shows the mean (scaled to the real er- 
rors) of the best-models obtained from Monte Carlo sim- 
ulations, when the errors are increased by 40% (see text). 
The dashed lines show the related 95% level, indicating 
an allowed range of Vdrc = 420 ± 40 km/s. Middle: the 
likelihood — In £gc of the Globular Cluster dataset. Bot- 
tom: the likelihood — ln£pjv of the Planetary Nebulae 
dataset. The models shown in Fig. |^ are circled. The 
dotted lines show the mean likelihood of the best-models 
obtained from Monte Carlo simulations. The dashed lines 
show the 95% level. 



our simple models based on logarithmic halo poten- 
tials fail to follow the knee inferred from the X-ray 
data, which Ikebe et al. (1996) interpret as the tran- 
sition of the galaxy potential to the Fornax cluster 
potential. Neverless these simple mass profiles are on 
the right mass scale even at these large distances. 

The integrated mass-to-light ratio in the B band 
is approximately constant at M/Lb ~ 10Mq/Lq out 
to i? w 60", increases to 12-15 at i? w 97" and reaches 
M/Lb = 22 - 48 Mq/Lq at i? = 417". Old {t > 8 
Gyrs) and metal rich (Z > Zq) stellar populations 
(Worthey 1994) can provide M/Lb ~ 10, matching 
at the same time the values of the Mg2 index and 
{B — V) color observed in the central regions of NGC 
1399 (Mg2 = 0.334, {B - V) = 0.99, Faber et al. 
1989). 

5. Conclusions 

We have presented new accurate and extended stel- 
lar kinematic data for NGC 1399, the cD galaxy of the 
Fornax cluster. Using a non-parametric algorithm, we 
have constructed smooth (at the 15 % level) spherical 
distribution functions which reproduce the observed 
velocity dispersion and /14 profiles at the four percent 
level. The models were tested for a grid of gravi- 
tational potentials, which combine a luminous part 
obtained from the deprojected surface brightness dis- 
tribution with the assumption of constant M/L, and 
a quasi-isothermal logarithmic dark halo potential. 
We have also explored the influence of a central black 
hole of mass 5 x IO^A/q , finding that the steep increase 
in (T in the inner 5 arcsec can be better followed. We 
have further constrained the mass profile of the galaxy 
by taking into account the likelihoods of the radial 
velocities measured for the galaxy's globular cluster 
and planetary nebulae systems and by comparing the 
results with the potential derived from X-ray mea- 
surements. Our main conclusions are summarized as 
follows: 

1. The stellar orbital structure of NGC 1399 is 
moderately radial out to i? ^ 60", with (3 = 
0.3 ± 0.1, and is not yet well-constrained at 
larger radii by the presently available data. 

2. The luminous component of the galaxy dom- 
inates the dynamics out to i? w 60", with 
M/Lb = IOMq/Lq. Old and metal rich stellar 
populations reproduce this high value, matching 
also the measured values of the central {B — V) 
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color and Mg2 index. 
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Fig. 7. — Top: the total mass as a function of radins for 
the range of acceptable models for NGC 1399. Bottom: 
the mass to light ratio in the B band. Line styles as in 
Fig. 1^ The full thick lines show the mass profiles derived 
from X-ray observations (see ^4.3). 



3. The influence of the dark halo is detected around 
i? w 60 — 90" and becomes dominant in the re- 
gion probed by the GC and PN systems, where 
the self-consistent model is ruled out. 

4. The circular velocity at the last measured stellar 
kinematic point is Vdac = 380 — 450 km/s. The 
corresponding mass profile range matches the 
mass curves derived from X-ray data with dif- 
ferent assumptions to model the galaxy's cool- 
ing flow. Inside R = 417" or 35 kpc we de- 
rive M = 1.2 — 2.5 X lO^^Af0 and a cumulative 
M/Lb = 22-48 Mq/Lq. 

5. The kinematics of the GCs and the PNs are con- 
sistent with those of the stars within the large 
statistical errors. The GCs appear slightly hot- 
ter, even when a correction for the different den- 
sity profile slope is taken into account. 

When compared to the normal giant ellipticals 
where a detailed stellar dynamical analysis has been 
performed (NGC 1600, Matthias and Gerhard 1999; 
NGC 2434, Rix et al. 1997 ; NGC 6703, G+98), 
the cD galaxy NGC 1399 appears to have a simi- 
lar anisotropy profile inside the inner few kiloparsecs, 
with (3 w 0.3. The luminous component of the galaxy 
dominates the gravitational potential out to 8 kpc 
from the center. As judged from the GC and PN 
data, the transition to the dark matter dominated 
regime appears to be more rapid than in the case of 
NGC 6703. It is possible that this difference reflects 
the special position of NGC 1399 at the center of the 
Fornax cluster. The best way to derive more stringent 
constraints on the nature of the transition region ap- 
pears to be from larger samples of PN velocities such 
as can be obtained with the VLT and other large tele- 
scopes. 

The present analysis based on spherical models is 
adequate to study the very round inner regions of 
NGC 1399 considered here. For more flattened sys- 
tems the methods developed by Matthias and Ger- 
hard (1999) and Cretton et al. (1999) are the future 
ways to follow. 
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